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We describe the relation between the isothermal atomic compressibility and density fluctuations 
in mixtures of two-component fermions with population or mass imbalance. We derive a generalized 
version of the fluctuation-dissipation theorem which is valid for both balanced and imbalanced Fermi- 
Fermi mixtures. Furthermore, we show that the compressibility, its critical exponents, and phase 
boundaries can be extracted via an analysis of the density fluctuations as a function of population 
imbalance, interaction parameter or temperature. Lastly, we demonstrate that in the presence of 
trapping potentials, the local compressibility and local density-density correlations can be extracted 
via a generalized fluctuation-dissipation theorem valid within the local density approximation. 
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Very recent experimental advances in Bose and Fermi 
systems have allowed for studies of density fluctuations 
and the use of the fluctuation-dissipation theorem to 
obtain information about some thermodynamic prop- 
erties of ultra-cold atoms. In the fermion case, the 
measurement of the density fluctuations and the atomic 
compressibility was extracted for non-interacting three- 
dimensional systems in harmonic traps [H while in 
the boson case, the connection between density fluctua- 
tions and compressibility was used to study superfluidity 
in a two-dimensional system, and extract critical expo- 
nents associated with the transition from a normal Bose 
gas to a Berezinskii-Kostcrlitz-Thouless supcrfluid 
The experimental extraction of the isothermal compress- 
ibility from density-fluctuation measurements were sug- 
gested several years ago both in harmonically confined 
systems [1] and optical lattices (Hj , but only recently im- 
provements in the detection schemes of density fluctua- 
tions became sufficiently sensitive to extract this infor- 
mation from experimental data [l|-Q . 

In principle, there are no major technical impediments 
to use the same technique for the study of density fluc- 
tuations in population imbalanced Fermi-Fermi mixtures 
of equal masses d, Q or unequal masses Q , where the 
compressibility and spin susceptibility matrix elements 
could be directly extracted from the density and density 
fluctuation profiles as discussed below. 

In this manuscript, we show that some local and global 
thermodynamic properties can be extracted through 
local measurements of densities and density fluctua- 
tions. We derive a generalized version of the fluctuation- 
dissipation theorem valid for any mixtures of atoms, and 
use it to analyze density fluctuations and the compress- 
ibility of mixtures of two-component fermions with and 
without population imbalance at low temperatures. For 
spatially uniform systems, we show that the global com- 
pressibility, phase boundaries and critical exponents can 
be extracted from measurements of the density and den- 
sity fluctuations for each component as a function of 



population imbalance, interaction parameter or temper- 
ature. While for spatially non-uniform systems, we show 
that the local and global compressibility can be extracted 
from measurements of the local density and local density 
fluctuations for each component. 

Hamiltonian: To investigate the physics described 
above, we start with the real space Hamiltonian {h = I) 
density for three dimensional s-wave supcrfluids 



H{r)^J2^i{r)(— fx^ + (r) ) (r) + ;7(r) , 
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(1) 

where ;7(r) = / dr'yi„t(r, r')V'|(r')V'l(r')V';(r)V't(r) con- 
tains the interaction potential T4nt(r, r') = —g5{r — r'), 
and creates fermions of mass nia labeled by index 

a. In addition, Va{v) and /i^ represent the trapping po- 
tential and the chemical potential for each fermion type, 
respectively. With the Hamiltonian H = J dr'H{r), we 
can study Fermi systems of equal masses = mi — m 
with population imbalance, or more generally we can 
study mixtures of fermions of unequal masses m-|- ^ toj,. 

From the grand partition function Z = Tre^^^^ , we 
can write the thermodynamic potential 57 = —TlnZ. 
First, we will ignore the trapping potential ^^(r) and 
discuss the spatially homogeneous case to simplify the 
discussion, but will return to the spatially inhomoge- 
neous situation later. We set ^^(r) = in Eq. ([1]) 
and rewrite the Hamiltonian H as Hi — X^aMa-^a, 

where = J dr V4(r) {~^) Mr) + U{r)] , 

and Na = J c?rV'i(r)^Q(r) is the number operator for 
hyperfine state a. The average number of particles 
Na — {Na) in hyperfine state a, defined by the ther- 



Nae~"/^ 



can be 



modynamic average (Na) 
rewritten in terms of the thermodynamic potential as 

Na = (Na) = -^n/^^la\T■ 

Pseudo-compressibility matrix: Next, we define the 
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pseudo-compressibility matrix as 



^af3 



T- 



T- 



(2) 



which through derivatives of Z can be expressed as the 
thermodynamic average = (TVaiV/j) — (7Vc<)(iV^). The 
mechanical stability of the system is guaranteed when 
both eigenvalues of Ka/3 are positive definite. Further- 
more, kap is a measure of density-density fluctuations: 



(3) 



The corresponding generalized compressibility matrix 
Kai3 can be obtained from kap through the relation = 

Kap/ {Na){Np) , describing a generalized fluctuation- 
dissipation theorem for multicomponent fermions. 

Using similar experimental techniques to those de- 
scribed in Refs. [H, it may be possible to measure 
the matrix elements of kap directly. Thus, it is impor- 
tant to identify the relation between the isothermal com- 
pressibility = —V{dP/dV)\T and the elements of the 
compressibility matrix Hap- 

Isothermal compressibility: The relation between kt 
and can be established by recalling that the ther- 
modynamic potential fl = —PV, where P is the pres- 
sure and V is the volume of the system. Defining 
G = ft + PV = 0, and recalling that $7 is a function 
of temperature T, volume V and chemical potentials /Xq, 
results in dG = -SdT - N■^d^if - N^d^i + VdP = 0. 
At constant temperature dT = 0, we can establish the 
relation —VdP\T = N^d^^lx + Nidiiilr- 

This means that the inverse isothermal compressibility 
= —V{dP/dV)T can be writen in terms of isother- 
mal partial derivatives of with respect to volume 
= N^{dii^/dV)\T + N^{dii^/dV)\T. But in turn the 
partial derivatives d^a/ dV\T can be expressed in terms 
dfia/dNplr and Np, leading to 
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T 
V 



TVyVt 



hi 



(4) 



This expression can be written in the compact form 



by using the definition Hap 



Hap/iNaNp). Therefore direct measurements of K^p lead 
to the isothermal compressibility kt of the system. 

Connection to pseudo-spin susceptibility: We can also 
work with the total number of particles = + N^, 
the particle number difference N_ = N-^ — N^_, and their 
corresponding chemical potentials ^± = {^^ ± /i4,)/2, 
respectively. In this case, we can define a similar 
pseudo-compressibility tensor kij = Td{Ni) /dfijlx — 
^T^^^l/^f^^^lJ\T, leading to S,, = {N,Nj) - (iV,)(7V,), 
where the indices i^j can each take ± values. The corre- 
sponding expression for kt has exactly the same form as 
before: 



1 



T 
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Nl 
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It is clear that kt reduces to the standard result 
for balanced populations where = and 7V+ = 
N, Kyi = TV-^Nl/k++, leading to the standard 
form of the fluctuation-dissipation theorem: kt = 



VT~ 



(iV2) - {Nf 



I (N)'^. The connection between the 
two representations is simple. The first diagonal term is 
K++ — K^f + k^i + 2k^^, the second diagonal term is 

K = k-f^ + ki^ ~ 2k-|-^, while the off-diagonal terms 

kjf = k |_ = K-f^i — ki^ are identical by symmetry. 

The pseudo-spin susceptibility Xzz of the system can 
also be extracted from the pseudo-compressibility matrix 
kij, since ^_ plays the role of an effective magnetic field 
hz along the quantization axis z, and A^- plays the role 
of the magnetization mz . For imbalanced Fermi systems 
this implies that 



= T- 



T dhz 



(6) 



and generalizes the results obtained for ultra-cold 
fermions interacting via p-wave interactions 0] . 

Compressibility matrix in a trap: In the presence of 
a trapping potential Va{r), we define the local chemi- 
cal potential as = A'a ^ Ki(r), and rewrite the 
Hamiltonian explicitly as H = Hi — J dr/iQ(r)nc(r), 
where ha{r) = ip^ {r)'ip{r) corresponds to the particle 
density operator. In this case, the grand partition func- 
tion is Z [T, V, Hair)] = Tre-t^i-/ drM„«„(r)]/T^ r^j^g ^^ler- 
modynamic potential [T, V, fia (r)] is also a functional 
of /ia(r), such that local and even non-local quanti- 
ties can be extracted. For instance the local density 
na(j) = (fiaij)) is simply written as the functional 
derivative ^^^(r) = Na{r) = —5^1/ 5 flair) , where Nair) 
is the local number of particles. Correspondingly the 
non-local pseudo-compressibility matrix is kapir,v') = 
—TS'^n/SiJ,air)SiJ,pir'), which in terms of the local parti- 
cle number operators N^ir) and iV^(r') becomes 

kapir,r') = {Nair)Npir')) - (iV„(r))(7V^(r')). (7) 

Since Nair) = (Nair)) is a local thermodynamic av- 
erage, we rewrite KQ^(r,r') = {6Nair)SNpir')) , where 
SNair) = Nair) — Nair) is the local fluctuation in 
the number of particles of type a. The local pseudo- 
compressibility matrix is simply kapir) — kapir,r). 

Within the local density approximation (LDA) the lo- 
cal isothermal compressibility kt ir) can be derived from 
the local pressure P(r) as Krir) — -~VdPir)/dV\T- Fol- 
lowing the steps leading to Eq. (HJ, we obtain 
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KT(r) 



1 



T 



Kapir)' 



(8) 



which is the local generalization of the fluctuation- 
dissipation theorem within LDA. Here, Kapir) = 
kapir)/ [Nair)Npir)]. 
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Imbalanced Fermi-Fermi mixtures: As an example of 
the general relations just derived we discuss the case 
of imbalanced Fermi- Fermi mixtures with equal masses, 
which has attracted a lot of interest lol- 15|. 

The resulting action corresponding to the Hamiltonian 
given in Eq. ^ has been succesfully calculated for the 
case of a uniform superfluid [l6| in the Gaussian ap- 
proximation as Sg ^ So + 5y A^('7)F~^A((7), where 
q = (q, I'f), with lyg ~ 2tt£T being the bosonic Matsub- 
ara frcqucny at temperature T. Here, A{q) is the order 
parameter fluctuation field and the matrix F~^{q) is the 
inverse fluctuation propagator. Furthermore, 



is the saddle point action, where i?k,Q — y + + l^kp + 
•SaCk,- is the quasiparticle energy when = f and is the 
negative of the quasihole energy when si = —1. We also 
use the notation £'k.± = (^k,t ^ ^k,4,) /2- In addition, 
Ak = AgFk is the order parameter for superfluidity for 
pairing with zero center of mass momentum, Fk = 1 for 
s-wave pairing, np{E]i^a) is the Fermi distribution, and 
Ck,± = (^k,t ± ^k,i) /2 = k^/{2m±) - fi±, where m± = 
2m-fm_i/ (mi ± TOf). Notice that m+ is twice the reduced 
mass of the t and ]. fermions, and that the equal mass 
case (m-t- = m^) corresponds to — > oo. 

The fluctuation term in the action leads to a cor- 
rection to the thermodynamic potential, which can be 
written as flc = fio + f^fluct, where fig — TSq and 
na^ct = TY^^lndet [TF-i(g)] . The saddle point condi- 
tion SVlo /(5 Aq = leads to the order parameter equation 
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"A'k.H 
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where Xk,± = (Ak,t ± -'^k,;) /2, with Xk,a = 
tanh (i?k.Q/(2T)) . As usual, we eliminate g in fa- 
vor of the scattering length via the relation 
l/.g = -m+y/(47ra,) + Ek |rk| V(2ek,+), where ek,± = 
(ek,t i ^k^) /2 = k"^ /{2m±). The order parameter equa- 
tion needs to be solved self-consistently with the number 
equations Na = ~7r^ ^ which has two contributions 



(10) 



where iVo^„ = -^t 
tion given by 



is the saddle point number equa- 



^0,a = 



1 - SaX\, 



2E] 



X 



k,+ 



kH 



(11) 



and iVfiuct,a = -9rifluct/9/^Q|T is the fluctu- 
ation contribution to Na given by A^fiuct a = 
-TEja [det¥-\q)\ / d,ia] / det¥-\q). 



To calculate Kq/j and kt from the thermodynamic 
potential Vt, we note that [Ao(/i-t', /ij.), //-|-, /xj^, T] — > 
VL /^|, T] . Because of the implicit dependence of Ag 
on and /ij^, the calculation of k^p requires 



= T 



T 



aiAol 



9|Aol 



(12) 



where the label "e" ("«") means explicit (implicit) deriva- 
tive. The mechanical stability of the uniform superfluid 
and normal phases is guaranteed if all eigenvalues of kap 
(or kij) are positive. This is achieved when Trit > and 
detK > 0. When the lowest eigenvalue of k reaches zero 
then the system becomes mechanically unstable. 



We define the Fermi momentum fc|n_|_ 



1,3 



l3 



where kpa is the Fermi momentum of each species, and 
the Fermi energy eF+ = fcF+/(2™+)- Foi' ^ mixture of 
fermions of equal masses, different hyperfine states and 
no trapping potential, the zero temperature phase dia- 
gram of population imbalance P — {N^ — N^)/{N^ + 
N^) = N-/Nj^ versus scattering parameter 1/(^^+0,,) 
is shown in Fig. [1^, where the normal (N), non-uniform 
(NU) and uniform (U) superfluid regions are indicated. 
The isothermal compressibility kt is shown in Fig. [TJd 
for \/{kFj^as) = 2.16 and changing P. Notice that as 
P increases, the dimensionlcss compressibility kt^f+Z^ 
diverges at a critical population imbalance Pc = 0.77, 
and becomes negative for P > Pc signaling a quantum 
phase transition from uniform superfluidity with coex- 
istence of excess unbound fermions and paired fermions 
in the same spatial region to a phase separated regime 
where excess unbound fermions and paired fermions tend 
to avoid being in the same region of space. 
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FIG. 1: a) The zero temperature phase diagram of popu- 
lation imbalance P = N-/N+ versus scattering parameter 
l/{kF+as) for equal mass fermions is shown, b) The di- 
mensionlcss isothermal compressibility kt^f^/V is shown for 
fixed l/{kF+as) = 2.16. 

In Fig. [2^, we show the dimensionlcss matrix ele- 
ments ka/seF+fT as a function of l/{kF^as) on the 
BEC side [l/{kF^as) > 0] for fixed population imbal- 
ance P = 0.5. Notice that kajs diverges as (A — Ac)~^, 
where A = l/(A:F+as), and Ac = 1.9 is the critical in- 
teraction parameter. As seen in Fig. ^jp, when pop- 
ulation imbalance is changed in the BEC regime, e.g. 
l/(^F+as) = 2.16; the dimensionlcss pseudo-spin sus- 
ceptibility Xzz^F+ = k eF+/T diverges at the phase 
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boundary between the uniform superfluid and the non- 
uniform phases as XzziF+ ^ {P — Pc)~^- The negative 
values of the matrix elements Kapep^ /T and Xzz^f+ just 
indicate the region of non-uniform superfluidity, i.e., the 
region where uniform superfluidity is not mechanically 
stable. Thus, the generalized fluctuation-dissipation the- 
orem described in Eq. ^ allows for the extraction of 
critical exponents of density-density and pseudospin- 
pseudospin correlations accross phase boundaries. 



the trap, leading to Xzz(r) 7^ 0. 
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FIG. 2: a) Pseudo-compressibility matrix elements 
kapeF+/T as a function of l/{kF+as) for population im- 
balance P = 0.5; b) Pseudo-spin susceptibility x^^^f+ = 
K__eF+/T as a function of P for interaction parameter 
l/{kF+as) = 2.16. 



In the case of a non-zero trapping potential Va{v), 
we make use of LDA and obtain the local thermody- 
namic potential r2(r) = fi (r), /X4,(r)] from the ther- 
modynamic potential in the absence of a trap O [^-|-, /.tj,] , 
via the substitution — > /^q (r) = Ha — Va{v). This im- 
plies that in Eqs. ([9]) and ([T0| the order parameter Ao 
and the number of particles iV^ become functions of po- 
sition r via the position dependent chemical potentials 
/ic((r). As a result we have Ao(r) = Aq \^i^{y), ^i{y)] and 
iV,(r) Ar„[Ait(r),Aii(r)]. 

We consider harmonic trapping potentials Va{v) ~ 
7Q,r^/2, where 7a = rriaUJa, with Wq, being the trapping 
frequencies of fermion of type a. For the equal mass case, 
we show in Fig. |3^, the particle number profiles Na{r) 
and the order parameter Ao(r) as a function of dimen- 
sionless position r/rj-Pi where r^F is the Thomas-Fermi 
radius defined through the condition epj^ — 7+r^^/2, 
where 7+ = 7-t- + 74,- These spatial profiles show that 
superfluidity coexists with excess unpaired fermions, but 
the majority of excess unpaired fermions are pushed away 
from the center of the trap. In Fig. [3}), we show the spa- 
tial dependence of Kai3{r), from which the local corre- 
lation functions (iVQ.(r)/V^(r)) and the pseudo-spin sus- 
ceptibility Xzzi'r) = K /T can also be easily extracted. 

Within LDA, KQ^(r) exhibit a discontinuous jump at the 
position Tc where Ao(rc) — 0. In the superfluid region, 
local particle fluctuations reveal the extra correlations 
brought in by full pairing such that Xzz(r) = 0. Out- 
side the superfluid region local particle fluctuations show 
a decrease in particle-particle correlations, as the excess 
unpaired fermions are pushed away from the center of 
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FIG. 3: a) Particle number profiles Na{r)/N+, and order 
parameter Ao(r)/eF+; b) Matrix elements kai3^F+/T as a 
function of v/ttf, for population imbalance P — 0.6 and 
interaction parameter l/(feF+as) ~ 3.0. 

Summary: We derived a generalized fluctuation- 
dissipation theorem for Fermi-Fermi mixtures, which was 
used to extract thermodynamic information (compress- 
ibility, spin-susceptibility, phase diagrams and critical ex- 
ponents) from density and density-fluctuation profiles of 
imbalanced mixtures of equal or unequal masses. We 
discussed continuum systems with and without trapping 
potentials. Using the local density approximation, we ob- 
tained expressions relating the local compressibility and 
local spin susceptibility to the local fluctuations in parti- 
cle numbers. Lastly, we applied our results to the case of 
population imbalanced Fermi systems of equal masses. 
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